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Abstract. - We study the detection of a BCS transition within a sample of Lithium-6 atoms 
confined in a harmonic trap. Using the local density approximation we calculate the pair 
correlation function in the normal and superfluid state at zero temperature. We show that the 
softening of the Fermi hole associated with a BCS transition leads to an observable increase in 
the intensity of off-resonant light scattered from the atomic cloud at small angles. 



The experimental realization of Bose-Einstein-Condensation in dilute Alkali vapors 
stimulated a variety of experimental and theoretical work in the field of quantum gases. 
Recently, attention was also drawn to the special properties of ultracold fermionic gases. In 
particular the absence of a s-wave scattering interaction in a spin polarized system would 
allow to realize an almost perfect example of a noninteracting Fermi gas, of which first 
indications were observed in an experiment with ^"^K atoms ^|. In the case of an attractive 
interaction between different hyperfine levels as occurs in ^Li |3|, the ground state will be a 
BCS superfluid. Stoof et al. Q first proposed to look for this phase transition in a mixture 
of two hyperfine states of ^Li, where due to the extraordinary large s-wave scattering length 
the transition temperature turns out to be in an observable range of about lOOnK. Their work 
was further expanded by Modawi and Leggett ||] , including all three hyperfine states that can 
simultaneously be enclosed in a magnetic trap. 

An experimental realization of the BCS-transition in magnetic traps is aggravated by par- 
ticular conditions regarding the life time of the sample, equal occupation of the hyperfine states 
and a large magnetic offset field to obtain the maximal scattering length of a = — 1140A for 
lithium so that cooling of the sample in a pure optical trap might be advantageous. Besides 
these difficulties, the most favorable way to detect a BCS-transition in future experiments is 
still a matter of discussion. In fact, the differences in the density and momentum distribution 
between the normal and superfluid phases are tiny so that unlike the case of bosonic gases, 
time-of-flight measurements will not be able to indicate the phase transition. As possible 
observable signatures collective excitations |^, changes in the line width and line shift in 
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light scattering experiments |10| and anomalous moments of inertia have been suggested. 
In this paper we show that the BCS transition in a cloud of cold atoms may be detected 
directly by off-resonant light scattering. This suggestion was made independently by Zhang et 
al. [p^ , however our results for the signature of the transition are quite different from theirs, 
see our discussion at the end of the paper. Recent measurements of the dynamic structure 
factor in Bose condensates show that even the time dependent pair correlation functions 
have indeed become accessible by experiments. 

We consider a system of fermionic atoms in two equally occupied hyperfine states trapped in 
an isotropic harmonic potential ^(x) ^muj'^x^ . Atoms in different hyperfine states interact 
via an attractive contact interaction ^(x, y) = Vo(5(x — y). The corresponding Hamiltonian is 



7? = ^ y d^x M/t (x)Ho*.(x) + Y E y rf'2;M/t (x)vl/t_^(x)*_.(x)*.(x). (1) 

where is a fermionic field operator which destroys a particle in hyperfine state a 

at location x and Tip = — ?i^/2mV^ + l^(x) — /i denotes the single particle Hamiltonian. 
The diagonalization of this Hamilton operator in mean field approximation is equivalent 
to solving the Bogoliubov-de Gennes equations, originally formulated for inhomogeneous 
superconductors Q 



Evu,,{^) = [Ho + VK(x)]?/k(x) + A(x)«k(x) (2) 
i^k«k(x) = -[?^o + VK(x)]«k(x) + A(xK(x). (3) 

Here we introduced the Hartree potential W{x) — Vo(^J-(x)\E'o-(x)) and the pair potential 
A(x) = — Vb(^'CT(x)^'_CT(x)), which both depend on the position in the trap. As in the 
conventional BCS theory the function vu{x) can be interpreted as the amplitude for the 
occupation of a Cooper pair with excitation energy i?k. The amplitudes are normalized by 



E K(x)<(y) + ^'k(x)«k(y)) - '5(x - y). (4) 

k 

We will not attempt to provide a complete solution of these equations which is quite 
involved numerically and also requires a nontrivial regularization procedure for the pseudo 
potential Instead, we use the local density approximation (LDA), which relies on the fact 
that both the Fermi wavelength A^;- and the BCS coherence length are much smaller than 
the typical scale xq = y^h/muj, on which the harmonic external potential varies. The results 
of the LDA can be obtained from the Bogoliubov-de Gennes equations by using an ansatz 
Mk(x) = exp(ikx)'Uk(x), Wk(x) — exp(ikx){jk(x) with slowly varying amplitudes WkCx), Wk(x). 
Inserting this into the Bogoliubov-de Gennes equations and neglecting all derivatives of 'Uk(x) 



and Wk(x)(^, the excitation energy is the familiar i!^k(x) — \/ (sk — ^^(x))^ + A2(x) with the 
local chemical potential /x(x) = /x — V{x) — W{x). Similarly the amplitudes take the standard 
BCS form 





(i)The LDA can also be interpreted as the semi-classical h ^ limit of an Ansatz tik(x) — 
^P''/''«k(x). 
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with the reduced single particle energies Ck(x) = Ek ^ order to determine the local 

values of /^(x) and A(x) we numerically solve the standard gap equation together with that 
for the total number of particles. The ultraviolet divergence in the gap equation which arises 
from the assumption of a contact potential with no extrinsic cutoff in energy of the attractive 
interaction, can be eliminated by using the standard relation between the bare interaction 
parameter Vq and the low energy effective interaction g = Anh'^a/m, which is determined by 
the scattering length a. The resulting density distribution and local pair potential are in very 
good agreement with those found earlier [Q, |ll|, |l5), with a typical cloud radius Rt ~ vp/uj. 
For particle numbers of order 6 • 10^ used in our calculations, the LDA is in fact an excellent 
approximation as was shown by Bruun et al. Jl5[ . 

We now turn to the calculation of the pair correlation function. Using the Wick-theorem 
the pair correlation function f;(x, y) separates into a normal part g^{x.,y) and an anomalous 
part g^(x,y): 



g(x,y) - (n(x)n(y)} - (n(x))(n(y)) - .g^(x, y) + .g^(x, y). (6) 

The normal part correlation function consists of the autocorrelation part and the anti- 
correlation between fermions in the same spin states. The anomalous correlation function is 
nonzero only in the superfluid phase, when anomalous expectation values occur due to the 
formation of pairs. At zero temperature and within the LDA the normal and anomalous parts 
of the correlation function are given by 

g^(x,y) = n(x)5(x-y)-2 (^y'^e«''(--'^)«k(x)^k(y)) ■ (7) 
<?^(x,y) = 2F(x,y)^^*(y,x). (8) 
where F(x, y) is the anomalous pair amplitude 



i^(x,y) = / ^ e'''(--'')^k(x)f)k(y). (9) 

The pair amplitude diverges for r = |x — y| ^ like mA(y)/47r?i^r, which is an artefact of 
the zero range interaction. Unlike the pair potential this divergence cannot be compensated 
by a divergence in the coupling constant. A rough estimate for the local pair correlator is 
obtained by taking F(x, x) — A(x.)/g with g as the renormalized interaction constant, leading 
to 



<?(x,x) = + ^. 10 

While the scattered intensity calculated below is rather insensitive to the precise value of 
(j((x, x), this relation shows explicitly that the positive anomalous correlations associated with a 
BCS-transition reduce the effect of the Fermi hole in the pair correlation function of fermions. 

In the normal state, where A(x) = 0, the anomalous correlations vanish and we obtain the 
following inhomogeneous generalization of a well-known result for Fermi gases; 



g-(x,y)=n(z)5(r)-- 



3 n{z) 
{kF{z)r) 



■(^sin{kF{z)r) — kp{z)r cos{kp{z)r) 



(11) 
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Fig. 1. - The scattering geometry used in tliis paper: A monodiromatic liglit beam witli wave vector k 
scatters at tiie atomic trap placed at the origin. In a minimal distance zq from the trap an observation 
screen is placed. The scattered intensity is measured at distance R± from the center of the observation 
screen. 

Here we introduced z — Max(a;, y), r = |x — y| and the local Fermi wave vector kp[z) — 
\J 2m^{z) / . The correlation function shows the typical Friedel oscillations with period Af/2 

and an algebraic decay r^'^ caused by the Fourier transform of the step in the Fermi distribution 
at zero temperature. In the superfluid regime with A(x) ^ 0, the pair correlations can only be 
calculated numerically. Reflecting the tiny differences of the density distributions, the normal 
correlation functions differ only at large distances, where their absolute values are negligible. 
The anomalous correlation function also exhibits Friedel oscillations and an exponential decay 
on a scale kp^. Naively one would expect a decay with the BCS coherence length = ^vp/^o 
as in the homogeneous case, but the inhomogenity of the amplitudes Ufc(x) and ^^(x) blocks 
this effect. Overall, the magnitude of the anomalous correlation function is about 10% of the 
normal correlations although it is formally a very small effect of order (Ao/ef)^ ~ 0.01 for 
realistic parameters in ^Li. 

In the following, we will see that due to an almost perfect cancellation of the autocorrelation 
and the normal correlation contributions in an off-resonance light scattering experiment, the 
effect of the anomalous correlations is quite appreciable. The scattering geometry is shown in 
Figure |l|. The incoming laser beam with amplitude E l and wave vector k produces an atomic 
polarization P(r) which for a large detuning S of the laser frequency from the resonance is 
given by 



Here d is the matrix element of the atomic dipole moment and n(r) the gas density in the 
trap. We calculate the scattered field E^c in first order Born approximation. Assuming that 
the incoming laser field is removed by the dark ground technique and neglecting the temporal 
variation due to the excitation processes, the measured intensity at position R± on the screen 



P(r) 



d2E(r) 
hS 



(12) 



is given by 1 16, 



/(R) = (E,,(R)E,,(R)} 



9/l 



(F„(i?^) + F,(i?^)), 



(13) 
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where is the laser intensity. The first contribution to the intensity is essentiahy the 
Fourier transform of the density distribution in the trap 



(14) 



where R and are unit vectors in the direction of R and in z-direction respectively. The 
second contribution to the intensity contains the pair correlation function: 



F,{R^)= / dVdVe*'=(^-®^>('— "'^g(r,r'). (15) 



The contribution Fn{R±) can be described as simple coherent scattering from a circular 
hole with an amplitude proportional to N'^ and indeed our numerical results show the typical 
Bessel function behaviour ||l^ . From the first minimum in the scattered intensity, which 
occurs at 1.22RX/d, we can extract an effective radius d/2 of the circular opening. Due 
to the variation of the density through the atomic cloud this radius turns out to be about 
half the Thomas-Fermi radius Rt of the trap. For the numerical calculations we use a trap 
frequency w = 27r x 144Hz and chemical potential = 110.5?iti; as typical parameters [Q. 
This corresponds to an oscillator length of xq ~ 3.4/im, a Thomas-Fermi radius of about 15a;o 
and a Fermi wavelength kp^ « 0.2x0 in the center of the trap. If we place the screen in a 
distance of zq = 2 cm from the trap, the first minimum in the scattered intensity is only 0.3 
mm away from the screen center. 

For the intensity contribution due to the correlations within the gas, the typical length 
scale is not the radius Rt of the trap, but the much smaller Fermi wavelength Xp. Since 
the scattering angle is proportional to the inverse typical length scale, the scattered intensity 
due to the correlations forms a much wider light cone. The Fermi wavelength is the typical 
scale for all parts of the correlation function. The normal part of the pair correlation function 
consists of two components, namely the autocorrelation function and the Friedel oscillating 
part which carry opposite signs. The positive autocorrelation leads to a quasi constant intensity 
proportional to N/R^, whereas the contribution of the Fermi hole with the Friedel oscillations 
is negative. The result is a strong suppression of the scattered intensity at small angles, 
reflecting the Fermi hole of the pair correlation function. While this effect could be used to 
detect the Fermi degeneracy in normal systems [^], the small positive contribution of the 
anomalous correlations also becomes significant. 

To demonstrate the effect of the BCS correlations, we numerically evaluated the total 
intensity on a screen with distance zq = 2cm from the center of the atomic cloud. The 
anomalous pair correlation function was determined by a refined Simpson routine, the integrals 



in (15), only one of which can be performed analytically, by a Monte Carlo method. We 
evaluated the pair correlation function at 200 million coordinate pairs (r,r') within the trap. 
The result is shown in figure |^, where the scattered intensity from a normal fluid and a 
superfluid gas of ^Li atoms are compared. For small distances R±^ from the center of the 
screen the superfluid phase transition significantly rises the intensity. The inset shows the 
'Fermi hole', the drastic reduction of the observed intensity of the fermionic gas in comparison 
with the intensity scattered from an ideal gas of atoms, where just the autocorrelation function 
contributes. The typical error bars on the curves as determined by the Monte Carlo method 
were below 5% for the normal part and below 10% for the anomalous part. 

As is shown in figure |3| the measured intensity of a normal fluid Fermi gas drops by almost 
90% in comparison with the ideal gas at zero temperature. This value for the normal fluid 
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Fig. 2. - The total scattered intensity distribution in arbitrary units with the distance R± in mm 
from the center of the observation screen. The dashed hne denotes the normal fluid state, the solid 
line the superfluid state. The inset also shows the intensity scattered from an ideal classical gas. 




Fig. 3. - The ratio of the measured intensities for different states with the distance R_i_ in mm from 
the center of the screen. The dashed line denotes the ratio of normal fluid fermionic gas to an ideal 
classical gas, the solid line the ratio of superfluid fermionic gas to the normal fluid gas. 



gas, however, is raised by almost 40% by the superfluid transition at small scattering angles, 
reflecting the softening of the Fermi hole by the anomalous correlations as in eq.(^0|). The 
peaks that occur within this graph result from side maxima of the central coherent scattering 
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peak, which for our system is about seven orders of magnitude stronger than the typical 
intensity caused by the correlations. Although the rise in the intensity will be smaller at finite 
temperatures, it should still be measurable. It is remarkable how the anomalous correlations 
qualitatively change the behaviour of the scattered intensity: at a fixed position one would 
first observe a drop in the intensity with falling temperature caused by Fermi statistics and 
below Tc a subsequent rise. 

Very recently, the suggestion to detect a BCS transition via off-resonant light scattering has 
also been made by Zhang et al. [|l2| . They found that the scattered light cone widens from an 
angle of order X/Rt in the normal state to an angle of order X/Xp in the presence of anomalous 
correlations. However in this work the normal part 5^(x, y) of the pair correlation function was 
neglected and the authors assumed R « R± to perform the integration over the ^-coordinate 
in equations (jlj) and (|l5|). Therefore the cancellation between the autocorrelation and the 
normal part of (;(x, y) which is responsible for the Fermi hole and its subsequent softening by 
the BCS correlations is not seen in their results. 

In conclusion we have performed a quantitative study of the possibility to detect a superfluid 
transition within a fermionic gas, e.g. ^Li, by measuring the intensity of off-resonant light 
scattered from the atomic trap. The intensity outside the central coherent scattering cone 
is dominated by the contributions of the pair correlation function of the system. The pair 
correlation function of a superfluid gas contains an anomalous part, which rises the observed 
intensity by as much as 40% compared to the degenerate Fermi gas, thus providing a clear 
signature for the onset of superfluid correlations. 
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